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Direct measurements of Na(p, 7) 3 Mg resonances and consequences for Na 

production in classical novae 
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The radionuclide 22 Na is a potential astronomical observable that is expected to be produced in 
classical novae in quantities that depend on the thermonuclear rate of the 22 Na(p, 7) 23 Mg reaction. 
We have measured the strengths of low-energy 22 Na(p, 7) 23 Mg resonances directly and absolutely 
using a radioactive 22 Na target. We find the strengths of resonances at E p = 213, 288, 454, and 610 
keV to be higher than previous measurements by factors of 2.4 to 3.2, and we exclude important 
contributions to the rate from proposed resonances at E p = 198, 209, and 232 keV. The 22 Na 
abundances expected in the ejecta of classical novae are reduced by a factor of ~ 2. 



PACS numbers: 29.30.Kv, 29.38.Gj, 26.30.Ca, 27.30.+t 



A classical nova is a thermonuclear outburst on the 
surface of a white-dwarf star that is accreting hydrogen- 
rich material from a binary companion. Such novae are 
ideal sites for the modeling of explosive nucleosynthe- 
sis because most of the relevant thermonuclear reaction 
rates are based on experimental information 0, Q ■ As- 
tronomical observations of novae enriched in neon at ul- 
traviolet, optical, and infrared wavelengths reveal that 
chemical elements as heavy as calcium can be synthe- 
sized and ejected Q, consistent with a model where the 
accreted layer mixes with neon seed material in the un- 
derlying white dwarf to synthesize heavier elements via 
a network of proton- induced reactions and j3 decays Q ■ 

Orbiting gamma-ray telescopes have the potential to 
detect lines that characterize the decays of particular iso- 
topes produced in novae, such as 26 Al and 22 Na. While 
novae may be one of many sources for the steady-state 
Galactic 26 Al abundance 0-0], its long half life (720,000 
yr) precludes its correlation with particular events. More 
detailed information on the underlying physical processes 
can be provided by the detection of Galactic 22 Na [5( , for 
which novae are expected to be the primary source. The 
2.6 yr half life of 22 Na is short enough that it will be 
localized, yet long enough that it will survive beyond the 
opaque explosive phase (and hence be visible) [8| , assum- 
ing it is not destroyed during the explosion. 

Recent modeling of 22 Na production in novae has fo- 
cused on oxygen-neon (ONe) white dwarfs @, 
which are more massive than common carbon-oxygen 
(CO) white dwarfs and provide a source of neon Mod- 
els of ONe nova predict 22 Na yields that are slightly lower 
than the upper limits set by the non-observation of the 
characteristic (1275-keV) 22 Na line in neon- type novae 
in the Galactic disk [ljj] ■ A key variable in the compari- 
son of models to observations is the destruction of 22 Na 
in novae, which depends crucially on the 22 Na(p, 7) 23 Mg 
reaction rate jl|, HHH 11 1 . 

The thermonuclear rate of the 22 Na(p, 7) 23 Mg reaction 



at peak nova temperatures (0.1 < T < 0.4 GK) is dom- 
inated by narrow, isolated resonances, and the currently 
accepted rate is based on two previous sets of direct mea- 
surements using radioactive 22 Na targets III, 14 1. The 
first set [l3[ yielded strengths for resonances with labora- 
tory proton energies E p > 290 keV. The second set 14 1 
identified a new resonance at 213 keV and provided a 
measurement of its strength relative to strengths from 



the first set [13j, with the conclusion that it dominates 
the rate. In 2004, measurements of the 12 C( 12 C,n7) 23 Mg 
reaction [llj provided an indirect basis for possible new 
resonances including one at E p = 198 keV, which the 
authors speculated could supersede the contributions of 
all other 22 Na(p, 7) 23 Mg resonances. 

We have made direct, absolute measurements of key 
22 Na(p, 7) 23 Mg resonance strengths for E p < 610 keV 
and searched for the proposed 198-keV resonance at 
the Center for Experimental Nuclear Physics and Astro- 
physics (CENPA) of the University of Washington using 
22 Na targets implanted at the radioactive ion beam facil- 
ity TRIUMF-ISAC. Unlike Ref. [H, we bombarded all 
target atom s by scanning the beam uniformly over the 
target area [15| and integrating the excitation function 
for each resonance. Our method is not very sensitive to 
target non-uniformity and stoichiomctry or to the evo- 
lution of the target distribution due to prolonged bom- 
bardment. More details will be given in Rcfs. liana . 

We produced 22 Na by bombarding thick SiC targets 
with a 40-/xA, 500-MeV proton beam from the TRIUMF 
cyclotron. A surface ionization source and a high resolu- 
tion mass separator were used to yield a 10-nA, 30-keV 
beam of 22 Na + . Three 300- /^Ci targets of 22 Na were fab- 
ricated by rastering the ion beam over a 5-mm diame- 
ter collimator and into an OFHC copper substrate. To 
prevent sputtering of 22 Na by the intense proton beam 
during the (p, 7) measurements, two targets were coated 
with 20 nm of Cr using vacuum evaporation 18]. The 
Cr-coated targets proved to be much more stable than 
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the bare targets and were used to acquire the data for all 
resonances, except the proposed one at E p = 232 keV. 

Proton beams of w 40 p,A were produced at CENPA by 
accelerating protons from an ion source at the terminal 
of a tandem Van de Graaff accelerator. The extracted 
beam was analyzed using a 90° magnet monitored with 
an NMR probe. The beam line was equipped with a mag- 
netic raster, a target chamber with a dual cold-shroud 
system, three electrically isolated collimators in series 
(including one on a sliding ladder), an electron suppres- 
sor, and a water-cooled target mount. Pressure in the 
chamber was held at w 1.5 x 10~ 7 torr. At each reso- 
nance energy the beam was tuned, with the raster off, 
through 3- and 1-mm diameter collimators on the sliding 
ladder backed by electrically isolated beam stops until 
transmissions of > 95% and > 50%, respectively, were 
achieved. After tuning, an 8-mm collimator open to the 
target was inserted, and the raster was used to irradiate 
the implanted target area uniformly during data acquisi- 
tion. The beam current on the target was integrated to 
measure the total number of incident protons, Nj,. 

Two detection systems were positioned at ±55° to the 
beam axis. Each system consisted of a high purity Ge 
(HPGe) crystal surrounded by Pb shielding and scintil- 
lators for cosmic-ray rejection. Pb plates 26-mm thick 
were positioned between the target and the HPGe detec- 
tors to reduce the 511-keV 7-ray detection rate. Data 
acquisition was triggered by signals from the HPGe de- 
tectors, and HPGe signals in coincidence with scintillator 
signals were rejected. Dead time was monitored using a 
pulser. The beam position on the target was inferred for 
each event by measuring and recording the vertical and 
horizontal magnetic fields produced by the raster. 

Measurements were made on known 22 Na(p, 7) 23 Mg 



resonances pj, |14|, which we find at E p = 213, 288, 454, 
and 610 keV, and on the proposed resonances at 198, 
209 11 1 , and 232 keV [3] . For each resonance, the beam 
energy was stepped over a range of « 25 keV, taking into 
account energy losses in the Cr and Cu (ss 4 keV, total), 
and the depth distribution of 22 Na in the target (of order 
10 keV-FWHM). The raster amplitude was scaled by the 
square root of the beam energy to maintain a constant 
illuminated area on the target. For diagnostic purposes, 
measurements were also made on a thick, extended 27 A1 
target, on a thick 5-mm diameter 27 Al disk embedded 
in OFHC copper, and on an ion-implanted 23 Na target 
similar to the 22 Na targets. 

If the beam is swept uniformly over the implanted area 
of the target, the resonance strength ujj may be deter- 
mined using the equation, 
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where Y is the 7-ray yield at a beam energy E, J Y dE 
is the integral over the excitation function, Ai a b is the 



proton reduced de Broglie wavelength, m is the proton 
mass, M is the target mass, Nt is the number of target 
atoms, pb = ^^/(Q/e) is the beam density normalized 
to the total number of incident beam particles, and ^j- 
is the areal density of the beam. 

The proton beam energy was calibrated using thick- 
target excitation functions for well-known 27 Al(p, 7) res- 
onances at E p = 326.6, 405.5, 504.9, and 506.4 keV ^ . 
The 7-ray energy spectra were calibrated using the 
7357.84- and 5088.05-keV [H| 7 rays from the 406-keV 
27 A\(p, 7) resonance and their first-escape peaks. Both 
of these calibrations were used to identify 22 Na(p, 7) 23 Mg 
resonances, and they provided two independent measure- 
ments of the resonance and excitation energies. The 
weighted averages of these values were adopted (Table HJ}. 

Detector efficiencies were determined for our setup, 
including the lead shielding, from absolute measure- 
ments at E 1 — 1332 keV using a 60 Co source cali- 
brated to 1.7% (99% C.L.), and relative measurements 
using 24 Na(/?7) decay branches O&v = 1369,2754 keV) 
and 27 Al(p, 7) reaction branches [2l], [22| at E p = 663 
keV {E 1 = 7575,10451 keV) and E p = 992 keV 
(E-y = 1779,4742, 10762 keV), together with detailed ef- 
ficiency simulations using the Monte Carlo code PENE- 
LOPE 23]. The photopeak detection efficiency ranged 
from 4.7 x 10" 4 to 2.1 x 10" 4 in the relevant energy range 
of 4.2 to 8.2 MeV, with an uncertainty of 6%, which in- 
cludes an estimated 3% uncertainty arising from possible 
22 Na(p, 7) 23 Mg anisotropy. 

The yield for each 7-ray branch at each proton energy 
was determined by integrating the photopeak and (in se- 
lected cases) the first-escape peak. The integrated beam 
current, the detection efficiency, a 55 to 65% correction 
for dead time, and a background subtraction were incor- 
porated into this integration. The dead-time correction 
was verified by measuring 27 Al(p,7) thick-target yields 
with and without a strong 22 Na source in the vicinity 
of the detectors. J Y dE was determined from the exci- 
tation function for each potential primary branch using 
energy windows of width ~ 25 and ~ 40 keV for the 
two HPGe detectors, integrating over proton energy, and 
summing over the branches (Fig. [T]). A (10 ± 5)% cor- 
rection was applied for the 213-keV resonance using the 
known shape of the excitation function because it was 
not completely covered by our data (Fig. [TJ (c)). 7 rays 
with energies below w 4 MeV could not be measured due 
to intense backgrounds from the radioactive targets. 

The initial value of Nt was determined for each 22 Na 
target by an in-situ measurement of target activity us- 
ing the 1275-keV line. 22 Na loss due to sputtering dur- 
ing proton bombardment [3] was monitored by measur- 
ing strong resonances repeatedly and measuring residual 
22 Na activity in the beam line after the removal of each 
target. The amount of potential loss for each measure- 
ment varied but was always < 12% for Cr-coated targets. 

The quantity pb was measured by comparing the yields 
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TABLE I: Energies and strengths for 22 Na(p, 7) 23 Mg reso- 
nances inves tig ated in the present work (Q- value = 7580.53 
± 0.79 keV HEIH). Present uncertainties and limits are at 
the 68% C.L. Finite strengths are the sum of partial strengths 
from measured branches, reported in detail elsewhere [l6lll7j]. 



Present 


Present 


Previous 


Present 


E p (keV) 


£ x (keV) 


uj'y (meV) 


CJ7 (meV) 


198° 


7769" 


< 4.0" 


< 0.51 e 


209" 


7780 a 


0.05 6 


< 0.40 / 


213.4 ± 1.4 


7784.6 ± 1.2 


1.8±0.7 C 


c 7 +i.6 


232° 


780 l a 


2.2 ± 1.0 d 


< 0.67 s 


288.0 ±1.1 


7856.1 ± 1.0 


15.8±3.4 C 


39 ±8 


454.2 ± 0.8 


8015.3 ± 0.8 


68 ± 20 c 


166 ± 22 


609.8 ± 0.8 


8163.9 ± 0.8 


235 ± 33 c 


591^° 3 



"Nominal value from Refs. [ill, llSj (not observed). 
^Indirect estimate from Ref. Hill . 
c Direct measurement from Refs. Il.'l llll . 
^Indirect estimate from Ref. [l9| . 

e Based on a search for the 5055-keV (58%) branch 11 
^Based on a search for the 5067-keV (66%) branch [Q 



(MeV) 



9 Based on searches 
(30.0%) branches [J 



for the 7350-keV (66.4%) and 7801-keV 



FIG. 1: Excitation functions (left column) and corresponding 
summed 7-ray spectra (right column). Panels (a) and (b) are 
for E p = 610 keV, panels (c) and (d) are for E p — 213 keV, 
and panels (e) and (f) are for E p = 198 keV. SE indicates 
a single-escape peak, and the arrow in (f) indicates where a 
potential 7 ray would be expected. 



from the 27 A\(p, 7) reaction using the thick, extended 
27 A1 target and the 5-mm disk 27 A1 target. The res- 
onances at 406 and 992 keV were used, yielding two 
values of pb in excellent agreement, whose average was 
adopted as the central value. To estimate a systematic 
uncertainty for pi,, this measurement was supplemented 
by tests in which the raster amplitudes and collimator 
diameters were varied and by using a Monte-Carlo sim- 
ulation that modeled the transport of the beam through 
the final components of the beam line. The result is 
Pb = 2.58 ± 0.25 cm -2 , where the uncertainty incorpo- 
rates potential non-uniformities in pb and target density. 

Using Eq. [TJ values of W7 were extracted (Table Q}. 
We find the strengths of the known 22 Na(p, 7) 23 Mg reso- 
nances at E p = 288, 454, and 610 keV to be higher than 



previous measurements [13j, [14| by a factor of 2.4 and 



that of the 213-keV resonance to be higher by a factor 
of 3.2. Our branches [l6| are roughly in agreement with 
Ref. [3] ■ We set upper limits on the strengths of the pro- 
posed resonances at 198, 209 keV [IH, and 232 keV [HI] 
by searching for known branches. To extract upper lim- 
its, the shape of the excitation function from a reference 
resonance (E p = 454 or 610 keV) was fit to the data 
of each proposed resonance, after stretching (set by the 
ratio of stopping powers), shifting, and normalizing the 
reference shape. The shift and normalization parameters 
were allowed to vary, and the upper limit was extracted 



to be the 68% quantile of the probability density function 
corresponding to the normalization. 

Our direct limit for the 198-keV resonance improves 
upon the indirect limit from Ref. llj by a factor of 8. 
We have set the first direct experimental limit for the 
proposed 11| resonance at 209 keV. Our direct limit for 
the 232-keV resonance is lower than the finite value of 
2.2 ± 1.0 meV reported previously [3], based on the (3- 
delayed proton decay of 23 Al. 

We tested our experimental method by making an ab- 
solute measurement of the 512-keV 23 Na(p, 7) resonance 
strength using the ion-implanted 23 Na target. We mea- 
sured the strength to be 79 ± 17 meV, consistent with 
the standard value of 91.3 ± 12.5 meV Q. 

Contributions to the thermonuclear 22 Na(p, 7) 23 Mg 
reaction rate were calculated using the presently mea- 
sured resonance energies and strengths in Table U (Fig. [2] 
(a)). Our upper limit for the 198-keV resonance and 
our strengths for the previously established resonances 
show that it is the 213-keV resonance that dominates the 
rate at peak novae temperatures ((Fig. [5] (b)). The 288- 
keV resonance begins to make a substantial contribution 
near the highest nova temperature. Due to the higher 
strengths for the 213- and 288-keV resonances, the to- 
tal reaction rate is significantly higher than the currently 
accepted rate ,[3 by roughly a factor of 3 (Fig.[2](c)). 

Since the 22 Na(p, 7) 23 Mg reaction is not expected to 
have a significant effect on the total energy generation 
in novae, we use the results of one-zone post-processing 
network calculations [H, 0] to test the effects of our rate 
on 22 Na production in nova models. We estimate that 
the amount of 22 Na produced is reduced by factors of 2 
to 3 compared to the previous rate 13J, [lj|, depending 
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FIG. 2: Panel (a) shows contributions to the thermonuclear 
22 Na(p, 7) 23 Mg reaction rate from individual resonances la- 
beled by E p in keV, based on the present measurements. 
Panel (b) shows the fractional contributions of selected reso- 
nances to the total present rate. Upper limits (68% C.L.) are 
denoted by dashed lines with arrows in (a) and (b). Panel (c) 
shows the ratios of reaction rates and their lo uncertainty lim- 
its to the central values of the reaction rate deduced from data 
in Refs. [H, (Tj on a log scale. The rates from Ref. [13, Q3] 
and the present work are represented by dashed lines in the 
vertical- and diagonal-hatched regions, respectively. Includ- 
ing the limit for the 198-keV resonance increases the present 
upper limit to the dot-dashed line. 



on the nova model used and the mass and composition 
of the underlying white dwarf. Treating the mixing of 
zones with the hydrodynamic code SHIVA 12811 indicates 
that these factors will be diluted by « 25% jSHl- 

In conclusion, we have measured low-energy 
22 Na(p, 7) 23 Mg resonance strengths using radioac- 
tive ion-implanted 22 Na targets and a technique that 
is more reliable than those employed in previous mea- 
surements LJ, |l4|. We find the strengths of the key 
resonances that destroy 22 Na in novae to be higher than 
the previous measurements 13, HH by factors of 2.4 to 
3.2. The contributions of proposed resonances at 198, 
209, and 232 keV to the destruction of 22 Na are found 
to be small, or negligible, at peak nova temperatures, 



and we have reduced uncertainties in the reaction rate. 
In summary, our measurements show that 22 Na will be 
destroyed much more efficiently in novae than previously 
thought, significantly reducing the prospects for the 
observation of 22 Na via the characteristic 1275-keV 
7~ray line using the spectrometer SPI onboard the 
currently-deployed INTEGRAL mission [30| . 
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